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Abstract 

This correspondence studies the downlink transmission in a multi-cell system, where multiple base stations (BSs) 
each with multiple antennas cooperatively design their respective transmit beamforming vectors to optimize the overall 
system performance. For simplicity, it is assumed that all mobile stations (MSs) are equipped with a single antenna 
each, and there is one active MS in each cell at one time. Accordingly, the system of interests can be modeled by a 
multiple-input single-output (MISO) interference channel (IC), termed as MISO-IC, with interference treated as noise. 
We propose a new method to characterize different rate-tuples for active MSs on the Pareto boundary of the achievable 
rate region for the MISO-IC, by exploring the relationship between the MISO-IC and the cognitive radio (CR) MISO 
channel. We show that each Pareto-boundary rate-tuple of the MISO-IC can be achieved in a decentralized manner 
when each of the BSs attains its own channel capacity subject to a certain set of interference-power constraints (also 
known as interference-temperature constraints in the CR system) at the other MS receivers. Furthermore, we show that 
this result leads to a new decentralized algorithm for implementing the multi-cell cooperative downlink beamforming. 

Index Terms 

Beamforming, cooperative multi-cell system, interference channel, multi-antenna, Pareto optimal, rate region. 

I. Introduction 

Conventional wireless mobile networks are designed with a cellular architecture, where base stations (BSs) from 
different cells control communications for their associated mobile stations (MSs) independently. The resulting inter- 
cell interference is treated as additive noise and minimized by applying a predesigned frequency reuse pattern such 
that the same frequency band is reused only by non-adjacent cells. Due to the rapidly growing demand for high-rate 
wireless multimedia applications, conventional cellular networks have been pushed towards their throughput limits. 
Consequently, many beyond-3G wireless technologies such as WiMAX and 3GPP UMTS Long Term Evolution (LTE) 
have relaxed the constraint on the frequency reuse such that the total frequency band becomes available for reuse by all 
cells. However, this factor-one frequency reuse pattern renders the overall network performance limited by the inter-cell 
interference; consequently, more sophisticated interference management techniques with multi-cell cooperation become 
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crucial. Among others, one effective method to cope with the inter-cell interference in the cellular network is via joint 

signal processing across different BSs. In this correspondence, we study a particular type of multi-BS cooperation 

for the downlink transmission, where we are interested in evaluating the benefit in terms of network throughput by 

cooperatively optimizing the transmit beamforming vectors for different BSs each with multiple antennas. Notice 

that the problem setup of our work is different from that for a fully cooperative multi-cell system considered in, 

£-g-> ffl-O, where a central processing unit is assumed with the global knowledge of all the required downlink 

channels and user messages to jointly design the transmitted signals for all BSs. In contrast, our work focuses on the 

decentralized implementation of the multi-cell cooperative downlink beamforming assuming only the local message 

and neighboring-channel knowledge at each BS, which is more practical than implementing the full baseband-level 

coordination. It is worth noting that decentralized multi-cell cooperative downlink beamforming has been studied in 

to minimize the total power consumption of all BSs to meet with MSs' individual signal-to-interference -plus-noise 

ratio (SINR) targets, based on the uplink-downlink beamforming duality. In this work, we provide a different design 

approach for rate-optimal strategies in decentralized multi-cell cooperative beamforming. 

For the purpose of exposition, in this work we consider a simplified scenario, where each MS is equipped with a 

single antenna, and at any given time there is only one active MS in each cell (over a particular frequency band). 

Accordingly, we can model the multi-cell cooperative downlink transmission system as a multiple-input single-output 

(MISO) Gaussian interference channel (IC), termed as MISO-IC. From an information-theoretic viewpoint, the capacity 

region of the Gaussian IC, which constitutes all the simultaneously achievable rates for all users, is still unknown in 

general (H, while significant progresses have recently been made on approaching this limit [9 ]. Capacity-approaching 

techniques for the Gaussian IC in general require certain signal-level encoding/decoding cooperations among the users, 

while a more pragmatic approach that leads to suboptimal achievable rates of the users is to allow only single-user 

encoding and decoding by treating the interference from other users as additive Gaussian noise. In this work, we 

adopt the latter approach to study the design of cooperative transmit beamforming for the MISO-IC. Particularly, we 

focus on the design criterion to achieve different rate-tuples for the users on the Pareto boundary of the achievable 

rate region for the MISO-IC. Due to the coupled signal structure, the achievable rate region for the MISO-IC with 

interference treated as noise is in general a non-convex set[] which renders the joint optimization of beamforming 

vectors to achieve different Pareto-boundary rate-tuples a challenging task. Note that this problem has been studied 

in ifTUl . where for the special two-user case, it was shown that the optimal transmit beamforming vector to achieve 

a Pareto-boundary rate -pair for the MISO-IC can be expressed as a linear combination of the zero-forcing (ZF) and 

'it is noted that the non-convex rate region is obtained without time-sharing (convex-hull operation) between different achievable rate-tuples. 
With time-sharing, the achievable rate region will become a convex set. 



maximum-ratio transmission (MRT) beamformers. The rate maximization for the IC with interference treated as noise 
has also been studied in the literature via various "pricing" algorithms (see, e.g., iTTTTl and references therein), while in 
general the price-based approach does not achieve the Pareto-optimal rates for the MISO-IC. In lfT2~l . the maximum 
sum-rate for the Gaussian IC with interference treated as noise is characterized in terms of degrees of freedom (DoF) 
over the interference-limited regime. 

In this correspondence, we develop a new parametrical characterization of the Pareto boundary for the MISO-IC in 
terms of the interference-power levels at all receivers caused by different transmitters, also known as the interference 
temperature (IT) levels in the newly emerging "cognitive radio (CR)" type of applications |fT3l . We show that each 
Pareto-boundary rate-tuple can be achieved in a decentralized manner when each of the users maximizes its own MISO 
channel capacity subject to a certain set of IT constraints at the other users' receivers, which is identical to the CR 
MISO channel transmit optimization problem studied in [14] and thus shares the same solution structure. We derive 
new closed-form solutions for the optimal transmit covariance matrices of all users to achieve an arbitrary rate-tuple 
on the Pareto boundary of the MISO-IC rate region, from which we see that the optimal transmit covariance matrices 
should all be rank-one (i.e., beamforming is optimal)o Furthermore, we derive the conditions that are necessary for 
any particular set of mutual IT constraints across all users to guarantee a Pareto-optimal rate-tuple for the MISO-IC. 
Based on these conditions, we propose a new decentralized algorithm for implementing the multi-cell cooperative 
downlink beamforming. For this algorithm, all different pairs of BSs independently search for their mutually desirable 
IT constraints (with those for the MSs associated with the other BSs fixed), under which their respective beamforming 
vectors are optimized to maximize the individual transmit rates. This algorithm improves the rates for the BSs in a 
pairwise manner until the transmit rates for all BSs converge with their mutual IT levels. 

Notation: I and denote the identity matrix and the all-zero matrix, respectively, with appropriate dimensions. For a 

square matrix S, Tr(S), \S\, and S 1 ^ 2 denote the trace, determinant, inverse, and square -root of S, respectively; 

and S ^ means that S is positive semi-definite ||T6l . Diag(a) denotes a diagonal matrix with the diagonal elements 

given by a. For a matrix M of arbitrary size, M H , M T , and Rank(iW) denote the Hermitian transpose, transpose, 

and rank of M, respectively. E[-] denotes the statistical expectation. The distribution of a circularly symmetric complex 

Gaussian (CSCG) random vector with the mean vector x and the covariance matrix 5] is denoted by CN(x, and 

~ stands for "distributed as". C mxn denotes the space ofmxn complex matrices. denotes the Euclidean norm 

of a complex vector (scalar) x. The log(-) function is with base 2 by default. 

2 We thank the anonymous reviewer who brought our attention to [ 151, in which the authors also showed the optimality of beamforming to 
achieve the Pareto-boundary rates for the Gaussian MISO-IC with interference treated as noise, via a different proof technique. 



II. System Model 

We consider the downlink transmission in a cellular network consisting of K cells, each having a multi-antenna 
BS to transmit independent messages to one active single-antenna MS. It is assumed that all BSs share the same 
narrow-band spectrum for downlink transmission. Accordingly, the system under consideration can be modeled by a 
K-user MISO-IC. It is assumed that the BS in the A;th cell, k = 1, . . . , K, is equipped with Mk transmitting antennas, 
Mk > 1. The discrete-time baseband signal received by the active MS in the kth cell is then given by 

K 

Vk = hk k x k + ^2 h fk x i + Zk W 

where Xk 6 C MfcXl denotes the transmitted signal from the kth BS; h kk 6 C lxMfe denotes the direct-link channel 
for the fcth MS, while h^ k G C lxMj denotes the cross-link channel from the jth BS to the kth MS, j / k; and z k 
denotes the receiver noise. It is assumed that ~ CA/"(0, a k ), V7c, and all z^'s are independent. 

We assume independent encoding across different BSs and thus x k s are independent over k. It is further assumed 
that a Gaussian codebook is used at each BS and Xk ~ CM(0, Sk), k = 1, . . . , K, where Sk = E[xkX k ] denotes the 
transmit covariance matrix for the kth BS, with Sk € C MkXMk and Sk h 0. Notice that the CSCG distribution has 
been assumed for all the transmitted signalsjj Furthermore, the interferences at all the receivers caused by different 
transmitters are treated as Gaussian noises. Thus, for a given set of transmit covariance matrices of all BSs, Si, ... , Sk, 
the achievable rate of the kth MS is expressed as 

Rk(Si, ...,S K ) = log (l + f H S f h k \ 2 ) • (2) 
V Z^khfkSjhjk + alJ 

Next, we define the achievable rate region for the MISO-IC to be the set of rate-tuples for all MSs that can be 
simultaneously achievable under a given set of transmit-power constraints for the BSs, denoted by Pi, . . . ,Pk- 

n± |J l(ri,...,r K ):0<r k <Rk(Si,...,S K ),k = l,...,K\. (3) 

{S k }:Tr(S k )<P k ,k=l,...,K 1 J 

The upper-right boundary of this region is called the Pareto boundary, since it consists of rate-tuples at which it is 
impossible to improve a particular user's rate, without simultaneously decreasing the rate of at least one of the other 
users. More precisely, the Pareto optimality of a rate-tuple is defined as follows [10]. 

Definition 2.1: A rate-tuple (n, . . . , rj() is Pareto optimal if there is no other rate-tuple (r[, . . . ,r' K ) with (r[, . . . ,r' K ) 
(ri, . . . , tk) and (r[,. . . , r' K ) 7^ (n, . . . , tk) (the inequality is component-wise). 

In this work, we consider the scenario where multiple BSs in the cellular network cooperatively design their transmit 

covariance matrices in order to minimize the effect of the inter-cell interference on the overall network throughput. In 

3 It is worth noting that in 1171 the authors point out that the CSCG distribution for the transmitted signals is in general non-optimal for the 
Gaussian IC with interference treated as noise, since it can be shown that the complex Gaussian but not circularly symmetric distribution can 
achieve larger rates than the symmetric distribution for some particular channel realizations. 
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Fig. 1. Achievable rate region and Pareto boundary for a two-user MISO Gaussian IC with interference treated as noise. 

particular, we are interested in the design criterion to achieve different Pareto-optimal rate-tuples for the corresponding 
MISO-IC defined as above. 

It is worth noting that in prior works on characterizing the Pareto boundary for the MISO-IC with interference 
treated as noise (see, e.g., [10] and references therein), it has been assumed (without proof) that the optimal transmit 
strategy for users to achieve any rate-tuple on the Pareto boundary is beaniforniing, i.e., Sk is a rank-one matrix 
for all /c's. Under this assumption, we can express Sk as Sk = w^w^k = 1,...,K, where € C Af *= xl denotes 
the beamforming vector for the kth user. Similarly as in the general case with Rank(Sfc) > 1, the achievable rates 
and rate region of the MISO-IC with transmit beamforming (BF) can be defined in terms of Wk's. However, it is 
not evident whether the BF case bears the same Pareto boundary as the general case with Rank(S'/ c ) > 1 for the 
MISO-IC. In this work, we will show that this is indeed the case (see Section Hill ). Accordingly, we can choose to 
use the rate and rate-region expressions in terms of either SVs or ii^'s to characterize the Pareto boundary of the 
MISO-IC, for the rest of this correspondence. 

In the following, we review some existing approaches to characterize the Pareto boundary for the MISO-IC with 
interference treated as noise. For the purpose of illustration, in Fig. \T\ we show the achievable rate region for a 
two-user MISO Gaussian IC with interference treated as noise (prior to any time-sharing of achievable rate -pairs), 
which is observed to be non-convex. A commonly adopted method to obtain the Pareto boundary for the MISO-IC 
is via solving a sequence of weighted sum-rate maximization (WSRMax) problems, each for a given set of user rate 



weights, Hk > 0, \/k, and given by 

K 

Max VVj. log (l + 7fc (wi,... ,w K )) 



k=l 

s.t. H^fell 2 <P k , k = l,...,K (4) 
where 7^(wi, . . . , iur-) is the receiver SINR for the kth user defined as 

7k ( W1 ,..., WK )= ^ fjri' — 5 , fc=i,... l ^. (5) 

It is easy to verify that this problem is non-convex, and thus cannot be solved efficiently. In addition, the above 
method cannot guarantee the finding of all Pareto-boundary points for the MISO-IC (cf. Fig. [TJ. 

An alternative method to characterize the complete Pareto boundary for the MISO-IC is based on the concept of 
rate profile lfT8l . Specifically, any rate-tuple on the Pareto boundary of the rate region can be obtained via solving the 
following optimization problem with a particular rate-profile vector, a = (ai, . . . , ax)'- 

Max. i? sum 

S.t. log(l+J k (wi,...,W K )) > OikRswn, k = l,...,K 

\\w k f<P k , k = l,...,K (6) 

with a k denoting the target ratio between the kth user's achievable rate and the users' sum-rate, R SU m- Without loss 
of generality, we assume that a k > 0, VA;, and Y2k=l a k = 1- F° r a given a, let the optimal solution of Problem © 
be denoted by Rg Um - Then, it follows that i?s Um ■ ct must be the corresponding Pareto-optimal rate-tuple, which can 
be geometrically viewed as (cf. Fig. Q]) the intersection between a ray in the direction of a and the Pareto boundary 
of the rate region. Thereby, with different q's, solving Problem ([6]) yields the complete Pareto boundary for the rate 
region, which does not need to be convex. 

Next, we show that Problem © can be solved via solving a sequence of feasibility problems each for a fixed r sum 
and expressed as 

Find {w k } 

s.t. log(l + 7 fe (xoi,...,xojr)) > a k r sum , k = l,...,K 

\\w k \\ 2 <P k , k = l,...,K. (7) 

If the above problem is feasible for a given sum-rate target, r sum , it follows that R* um > r sum ; otherwise, R* um < r sum . 
Thus, by solving Problem (O with different r sum 's and applying the simple bisection method |fl6l over r sum , i?* um 
can be obtained for Problem ©. Let */ k = 2 aferaum — 1, k = 1, ... K. Then, for Problem ([7]), we can replace the rate 
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constraints by the equivalent SINR constraints given by 

j k (wi,... ,w K ) > 7 fc , k = l,...,K (8) 

Similarly as shown in [19], the resultant feasibility problem can be transformed into a second-order cone programming 
(SOCP) problem, which is convex and can be solved efficiently [20]. 

III. Characterizing Pareto Boundary for MISO-IC via Interference Temperature Control 

In this section, instead of investigating centralized approaches, we present a new method to characterize the Pareto 
boundary for the MISO-IC in a distributed fashion, by exploring its relationship with the CR MISO channel fl4l . 
We start with the general-rank transmit covariance matrices Sk's for the MISO-IC. First, we introduce a set of 
auxiliary variables, Tkj , k = 1, . . . , K , j = 1, . . . , K, j ^ k, where T^j is called the interference -power or interference- 
temperature (IT) constraint from the feth BS to jth MS, j ^ k, T^j > 0. For notational convenience, let T be the 
vector consisting of all K(K — 1) different T^'s, and be the vector consisting of all 2(K — 1) different T^'s and 
F jk \ j = 1, . . . ,K,j ^ k, for any given k G {1, . . . , K}. 

Next, we consider a set of parallel transmit covariance optimization problems, each for one of the K BSs in the 
MISO-IC expressed as 

Max. log 1+ — j 

s.t. h^Skhkj < F h j, Vj / k 

Tr(-S fe ) < P k , S k y (9) 

where k G {1, . . . , K}. Note that in the above problem for a given k, is fixed. For notational convenience, we 
denote the optimal value of this problem as C^fT^). If in the objective function of (© we set T^ = h^Sjhj^^j ^ k, 
Cfc(rfc) becomes equal to the maximum achievable rate of an equivalent MISO CR channel lfl4l . where the Mi BS 
is the so-called "secondary" user transmitter, and all the other K — 1 BSs, indexed by j = 1,... ,K, j / k, are 
the "primary" user transmitters, each of which has a transmit covariance matrix, Sj, and its intended "primary" user 
receiver is protected by the secondary user via the IT constraint: h^Skh^j < Tf,j. In 031, it was proved that the 
solution for Problem (© is rank-one, i.e., beamforming is optimal, and in the special case of K = 2 (i.e., one single 
primary user), a closed-form solution for the optimal beamforming vector was derived. In the following proposition, 
we provide a new closed-form solution for Problem (© with arbitrary values of K, from which we can easily infer 
that beamforming is indeed optimal. 



Proposition 3.1: The optimal solution of Problem © is rank-one, i.e., S k = w k w¥ , and 



k 

w k= \Y1 X kjh k jh^ + A fcfc i" h kk y/Pk (10) 

where X k j, j 7^ fc, and A^fc, are non-negative constants (dual variables) corresponding to the kth BS's IT constraint 

for the jth MS and its own transmit-power constraint, respectively, which are obtained as the optimal solutions for 

the dual variables in the dual problem of Problem ©; and p k is given by 

( 1 Ej^jk + cri y 1 
h yln2 \\A k hkk\\ 2 J \\A k h kk \\ 2 

where A k = (y2j^ k X k jh k jh^ + X kk I^ 1 and (x) + = max(0, x). 

Proof: Please see Appendix [A] ■ 
Now, we are ready to present a parametrical characterization of the Pareto boundary for the MISO-IC in terms of 
r as follows. 

Proposition 3.2: For any rate-tuple . . . ,Rk) on the Pareto boundary of the MISO-IC rate region defined in 
(0), which is achievable with a set of transmit covariance matrices, Si, ... , Sk, there is a corresponding interference- 
power constraint vector, T > 0, with T k j = h^jS k h k j,\/j / k,j G {1, ...,K}, and k £ {1, . . . , K}, such that 
Rk = C k (T k ),\lk, and S k is the optimal solution of Problem © for the given k. 

Proof: Please see Appendix iBl ■ 

From Proposition 13.21 it follows that the Pareto boundary for the MISO-IC is parameterized in terms of a lower- 
dimensional manifold over the non-negative real vector T, in comparison with that over the complex transmit covariance 
matrices, S k s, or with that over the complex beamforming vectors, w k s. Furthermore, by combining Propositions 
13.11 and 13.21 it follows that beamforming is indeed optimal to achieve any rate-tuple on the Pareto boundary for the 
MISO-IC. 

Remark 1: It is worth noting that the dimensionality reduction approach proposed in this work for characterizing 
the Pareto boundary of the MISO-IC is in spirit similar to that proposed in [10], where it has been shown that the 
transmit beamforming vectors to achieve any Pareto-boundary rate-tuple of the K-user MISO-IC with interference 
treated as noise can be expressed in the following forms: 

K 

w k = ^2Ckjh k j, k = l,...,K (12) 

where £, k /s are complex coefficients. Note that under the assumption of independent h k fs, the above beamforming 
structure is non-trivial only when M k > K. For this case, from Remark [2] in Appendix |A] it is known that for the 
optimal beamforming structure given in (fTOl . we have \ kk > 0. With this and by applying the matrix inversion lemma 
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11211 . it can be shown (the detailed proof is omitted here for brevity) that the optimal beamforming structure given 
by (fTOb is indeed in accordance with that given by ([T2l . The main difference for these two methods to characterize 
the Pareto boundary for the MISO-IC lies in their adopted parameters: The method in our work uses K(K — 1) real 
Tfcj's, while that in iflOl uses K(K — 1) complex £fcj's. Note that Tkj corresponds to the IT constraint from the A;th 
user transmitter to the jth user receiver, whereas there is no practical meaning associated with Consequently, 
as will be shown next, the proposed method in our work leads to a practical algorithm to implement the multi-cell 
cooperative downlink beamforming, via iteratively searching for mutually desirable IT constraints between different 
pairs of BSs. 

IV. Decentralized Algorithm for Multi-Cell Cooperative Beamforming 

In this section, we develop a new decentralized algorithm that practically implements the multi-cell cooperative 
downlink beamforming based on the results derived in the previous section. It is assumed that each BS in the cellular 
network has the perfect knowledge of the channels from it to all MSs. Furthermore, it is assumed that all BSs operate 
according to the same protocol described as follows. Initially, a set of prescribed IT constraints in T are set over the 
whole network, and the A;th BS is informed of its corresponding Tfc, k = 1, . . . , K. Accordingly, each BS sets its own 
transmit beamforming vector via solving Problem © and sets its transmit rate equal to the optimal objective value 
of Problem ©, which is achievable for its MS since the actual IT levels from the other BSs must be below their 
prescribed constraints. Then, by assuming that there is an error-free link between each pair of BSs, all different pairs 
of BSs start to communicate with each other for updating their mutual IT constraints (the details are given later in this 
section), under which each pair of BSs reset their respective beamforming vectors via solving Problem (O such that 
the achievable rates for their MSs both get improved. Notice that each pair of updating BSs keeps the IT constraints 
for the MSs associated with the other BSs excluding this pair fixed; and as a result, the transmit rates for all the other 
MSs are not affected. Therefore, the above algorithm can be implemented in a pairwise decentralized manner across 
the BSs, while it converges when there are no incentives for all different pairs of BSs to further update their mutual 
IT constraints. 

Next, we focus on the key issue on how to update the mutual IT constraints for each particular pair of BSs to 
guarantee the rate increase for both of their MSs. To resolve this problem, in the following proposition, we first 
provide the necessary conditions for any given T > (component-wise) to correspond to a Pareto-optimal rate-tuple 
for the MISO-IC, which will lead to a simple rule for updating the mutual IT constraints between different pairs of 
BSs. Note that from Proposition 13.21 it follows that for any Pareto-optimal rate-tuple of the MISO-IC, there must 
exist a T such that the optimal solutions of the problems given in (O for all /c's are the same as those for the general 



to 

formulation of MISO-IC to achieve this rate-tuple. However, for any given T > 0, it remains unknown whether this 
value of r will correspond to a Pareto-optimal rate-tuple for the MISO-IC. 

Proposition 4.1: For an arbitrarily chosen T > 0, if the optimal rate values of the problems in © for all fc's, 
Ck(Tk)'s, are Pareto-optimal on the boundary of the MISO-IC rate region defined in ©, then for any pair of G 



{1, . . . , K},j € {1, ... , K}, and i ^ j, it must hold that |.D«| = 0, where D; t j is defined as 

(13) 



Proof: Please see Appendix O ■ 
Note that Difs for all different pairs of can be obtained from the (primal and dual) solutions of the problems 
given in © for all fc's with the given T (for the details, please refer to Appendix lAl. More specifically, we have 

dd {Ti] 



Ai 7 - (14) 



where \j is the solution for the dual problem of Problem (O with k = i, which corresponds to the jth IT constraint, 
and from the objective function of Problem ©, 



ar* ~ in2(E^r a + a?)(E^r a + a? + hf5*h i ,. 



(15) 



dC (T ) dC (r ■) 

where S* is the optimal solution of Problem (© with k = i. Similarly, — gr . an d — or can ^ e obtained from 
solving Problem (O via the Lagrange duality method with k = j. 

From Proposition 14- 1 1 the following observations can be easily obtained (the proofs are omitted for brevity): 

• For any particular T that corresponds to a Pareto-optimal rate-tuple, it must hold that Tij < f y , ^ j, 

\\h H -h-\\ 2 P- 

where Tij = — " — - corresponds to the case of using maximum transmit power with MRT beamforming for 



the ith BS; 

Hj ""IJ ~ 1 IJ> 



For any particular T that corresponds to a Pareto-optimal rate-tuple, it must hold that h^S*hij = T^, Vi, j, i ^ j, 



i.e., the IT constraints across all BSs must be tight. 

From the above observations, we see that if we are only interested in the values of T that correspond to Pareto- 
optimal rate-tuples for the MISO-IC, it is sufficient for us to focus on the subset of T within the set T > 0, in which 
Tij < Tij and r„ = h%S*hij,Vi,j,i^ j. 

Based on Proposition 14.11 we can develop a simple rule for different pairs of BSs in the cooperative multi-cell 
system to update their mutual IT constraints for improving both of their transmit rates, while keeping those of the 
other BSs unchanged. From the proof of Proposition 14.11 given in Appendix |Cj it follows that the method for any 
updating BS pair to fulfill the above requirements is via changing T^ and Tji according to (I38T ). Note that in 
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Fig. 2. Achievable rates for the proposed algorithm in a two-user MISO Gaussian IC with interference treated as noise. 

general, the choice for dij in ( |38T ) to make Dijdij > is not unique. For notational conciseness, let = 
it can then be shown that one particular choice for d^ is 

d^ = sign (ad — be) ■ [aijd — b,a — ctijc] T (16) 

where sign(x) = 1 if x > and = — 1 otherwise; otij > is a constant that controls the ratio between the rate 
increments for the ith and jth BSs. It can be easily verified that when > 1, a larger rate increment is resulted for 
the ith BS than that for the jth BS, and vice versa when aij < 1 (provided that the step-size 5ij in (I38T ) is sufficiently 
small). 

More specifically, the procedure for any BS pair (i,j), i ^ j, i € {1, . . . , K}, and j G {1, . . . , K}, to update their 
mutual IT constraints is given as follows. First, the 2th BS computes the elements a and b in according to (fT4l ) 
and ( fT5l ), respectively, with the present value of IV Similarly, the jth BS computes c and d with the present value 
of Tj. Next, the ith BS sends a and b to the jth BS, while the jth BS sends c and d to the ith BS. Then, assuming 
that and 5ij are preassigned values known to these two BSs, they can both compute d^ according to (fT6l ) and 
update Tij and Tji according to d38l ) in Appendix ICl Last, with the updated values V- and Y'^, these two BSs reset 
their respective beamforming vectors and transmit rates via solving (O independently. Note that the above operation 
requires only local information (scalar) exchanges between different pairs of BSs, and thus can be implemented at a 
very low cost in a cellular system. One version of the decentralized algorithm for cooperative downlink beamforming 
in a multi-cell system is described in Table U Since in each iteration of the algorithm the achievable rates for the 
pair of updating BSs both improve and those for all other BSs are unaffected (non-decreasing), and the maximum 
achievable rates for all BSs are bounded by finite Pareto-optimal values, the convergence of this algorithm is ensured. 




a b 
c d 
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Initialize T > in the network 

BS k sets Wk via solving $9$ with the given I\, k — 1, . . . , K 
Repeat 

For i = 1, . . . , K , j = 1, . . . , K, j i, 
BS i computes a and b in Dij (cf. d 14b . d!5l l) with the given Ti 
Similarly, BS j computes d and c in with the given Tj 
BS i sends a and 6 to BS j 
BS j sends c and d to BS i 

BS i/j computes dij (cf. l!16t). then updates I\y and r# (cf. ( 1381 )) 
BS i/j resets Wi/wj via solving with the updated Ti/Tj 
End For 
Until |D y | = 0,Vi / j. 

TABLE I 

Algorithm for Cooperative Downlink Beamforming. 



Example 4.1: In Fig. [2] (with the same two-user MISO-IC as for Fig. [[]), we show the Pareto boundary for an 
example MISO-IC with K = 2, Mi = M 2 = 3, Pi = 5, P 2 = 1. and a\ = a\ = 1, which is obtained by the 
proposed method in this correspondence, i.e., solving the problems given in (© for k = 1, 2, and a given pair of 
values Ti2 and I^i with < Ti2 < F12 and < T2i < T21, and then taking a closure operation over the resultant 
rate -pairs with all different values of Ti2 and I^i within their respective ranges. We demonstrate the effectiveness 
of the proposed decentralized algorithm for implementing the multi-cell cooperative downlink beamforming with two 
initial rate-pairs, indicated by "ZF" and "MRT" in Fig. |2j which are obtained when both BSs adopt the ZF and the 
MRT beamforming vectors, respectively, with their maximum transmit powers. It is observed that the achievable rates 
for both MSs increase with iterations and finally converge to a Pareto-optimal rate-pairo By comparing the two cases 
with Q12 = 1 and ai2 = 10, it is observed that a larger value of 012 results in larger rate values for the first MS in 
the converged rate-pairs, which is in accordance with our previous discussion. 

V. Concluding Remarks 

In this correspondence, based on the concept of interference temperature (IT) and under a cellular downlink setup, 
we have developed a new method to characterize the complete Pareto boundary of the achievable rate region for the 
iT-user Gaussian MISO-IC with interference treated as noise. It is shown that the proposed method also leads to a new 
decentralized algorithm for implementing the downlink beamforming in a cooperative multi-cell system to achieve 
maximal rates with a prescribed fairness guarantee. 

4 We have verified with a large number of random channels and a variety of system parameters that the proposed algorithm always converges 
to Pareto-optimal rate-pairs for the two-user MISO-IC with randomly selected initial rate-pairs. However, we could not prove this result in 



general by, e.g., showing that the conditions given in Proposition 14. 1 1 are not only necessary (as proved in this work) but also sufficient for any 
given r to achieve a Pareto-optimal rate-tuple for the MISO-IC. 



13 

There are a number of directions along which the developed results in this work can be further investigated. First, 
it would be interesting to extend the multi-cell cooperative beamforming design based on the principle of IT to the 
scenario where each BS supports simultaneous transmissions to multiple active MSs each with a single antenna or 
multiple antennas. Second, it remains yet to be proved whether the necessary conditions derived in this work for 
any particular set of IT constraints across the BSs to guarantee a Pareto-optimal rate-tuple for the MISO-IC are 
also sufficient, even for the special two-user case. This proof is essential for the proposed downlink beamforming 
algorithm to achieve the global convergence (Pareto-optimal rates). Last but not least, it is pertinent to analyze the 
proposed decentralized algorithm that iteratively updates the mutual IT constraints between different pairs of BSs 
from a game -theoretical viewpoint. 

Appendix A 
Proof of Proposition 13. II 

It can be verified that Problem © is convex, and thus it can be solved by the standard Lagrange duality method 
[TU)1. Let Afcj, j ytz k, and X kk be the non-negative dual variables for Problem (O associated with the kth BS's IT 
constraint for the jth MS and its own transmit-power constraint, respectively. The Lagrangian function for this problem 
can be written as 

where X k = [A^i, . . . , X k K\- The dual function of Problem (© is given by 



L(S k , A fc ) = log ( 1 + h ^ k h kk _ \ _ J2 X kj (h%S k h kj - T kj ) - \ kk (Tr(S k ) - P k ) (17) 



Accordingly, the dual problem is defined as 



g(Afc) = max L(S k ,\ k ). (18) 
S k yo 



min g(X k ) (19) 
A fc >o 



where A& > means component-wise non-negative. Since Problem Q is convex with strictly feasible points lfl6l . 
the duality gap between its optimal value and that of the dual problem is zero; thus, Problem ([9]) can be equivalently 
solved via solving its dual problem. In order to solve the dual problem, we need to obtain the dual function g(X k ) for 
any given X k > 0. This can be done by solving the maximization problem given in (fT8T ). which, according to (fTTT ). 
can be explicitly written as (by discarding irrelevant constant terms): 

Max. log ( 1 + J l * k r khk * g | " Tr(B k (X k )S k ) 

s.t. s k y o (20) 
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where B k (X k ) — J2j^ k ^kjh k jh k j + AfcfcJ and B k (X k ) ^ of dimension M k x M k . In order for Problem (l20l to 
have a bounded objective value, it is shown as follows that B k (Xk) should be a full-rank matrix. Suppose that B k (X k ) 
is rank-deficient, such that we could define S k = q k v k v^ > wnere Qk > and v k € C MfcXl satisfying \v^\ = 1 and 
B k (X k )v k = 0. Thereby, the objective function of Problem (l20l reduces to 

Hflt .?% V (21) 

Due to the independence of h kk and h k j's, and thus the independence of h kk and t?j., it follows that > 

with probability one such that (T2TT) goes to infinity by letting q k — > oo. Since the optimal value of Problem (O must 
be bounded, without loss of generality, we only need to consider the subset of Xk in the set A^ > to make B k {X k ) 
full-rank. 

Remark 2: Note that from the definition of B k {X k ) and the Karush-Kuhn-Tucker (KKT) optimality conditions |[T6l 
of Problem (©, it follows that B k {X k ) is full-rank only when either of the following two cases occurs: 

• Xkk > 0: In this case, the transmit power constraint for the kth BS is tight for Problem (|9]); 

• Afcfc = 0, but there are at least M k A^ 's, j ^ k, having Ajy > 0: In this case, regardless of whether the transmit 
power constraint for the kth BS is tight, there are at least M k out of the K — 1 IT constraints of the kth BS are 
tight in Problem (©. Note that this case can be true only when M k < K — 1. 

From the above discussions, it is known that (B k (X k ))^ 1 exists. Thus, we can introduce a new variable Sk for 
Problem (J20]) as 

S k = (BkiX^SkiBkiXk)) 1 / 2 (22) 



and substituting it into (1201 ) yields 

Max. log f 1 + <(B t (A t »-/^(ft(A t ))-V^A _ g 

s.t. S k t 0. (23) 

Without loss of generality, we can express Sk into its eigenvalue decomposition (EVD) as Sk = U k & k U^ , where 
U k = [u k i, ■■■ , u kMk ] £ c MkXMk is unitary and @ k = Diag([6» fcl , . . . , 9 kMk \) h 0. Substituting the ED of S k into 
(|23l yields 

i og f i + r 1/a « M ii a \ _ ^ 9h 

U k ,® k \ l^j^ k L jk + a k J i=1 

s.t. || ttfci || = ljVi, UfaUki = 0,\/l ^ i 

O kl >0,Vi. (24) 
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For any given U k , it can be verified that the optimal solution of 0^ for Problem d24b is given by 

- ,,^ ( 1: ( ao)^u M iiO + if * = argmax ^^'-'^> ii^c^(^))- i/a ««ii (25) 

otherwise. 
Thus, it follows that for the optimal solution of Problem (T2"3l . Rank(Sfc) < 1. Furthermore, let %' denote the index 

of % for which 6^1 > 0. The objective function of Problem (T2"4l ) reduces to 

/ e ki ,\\h¥,(B k (\ k ))- l i 2 u ki ,\\ 2 \ „ 

log 1+ k ^ fc V 2 fa " - (26) 

Clearly, the above function is maximized with any 9 k i> > when 

(B k (\ k ))-^h kk 
~ \\(B k (\k))-V 2 hkk\\' K ° 



From (1251) and (1271) . it follows that the optimal solution for Problem (1231 ) is 

„ ln ]J h ;: ( TyT7 (B k (x k ))-^h kk h^B k (x k) r^. (28 ) 

Combining the above solution and d22l) . it can be easily shown that the optimal solution S k for Problem (© is as 
given by Proposition 13.11 

With the obtained dual function g(X k ) for any given A&, the dual problem ( fT9l can be solved by searching over 
Afc > to minimize g(X k ). This can be done via, e.g., the ellipsoid method 11221 . by utilizing the subgradient of 
g(X k ) that is obtained from (fTTl) as T k j — h k jS* k h k j for X k j,j / k and P k — Ti(S* k ) for \ kk , where is the optimal 
solution for Problem d20l ) with the given A&. When Xf, converges to the optimal solution for the dual problem, the 
corresponding S* k becomes the optimal solution for Problem (©. Proposition 13.11 thus follows. 

Appendix B 
Proof of Proposition [32] 

Since the given set of Si, ... , Sk achieves the Pareto-optimal rate-tuple {R\, . . . , Rr) for the MISO-IC, from (f2]) 
and © it follows that for any k G {1, . . . , K } 

i4 = log(l + h *fkh kk A (29) 



Since Tj k = h^ k Sjhj k ,\/j ^ k, (1291 can be rewritten as 



h kk S k h kk 



fl fe = log 1 ™ • (30) 

Note that (1301 is the same as the objective function of Problem © for the given k. Furthermore, from ((3]) it follows 
that Tr(S k ) < P k . Using this and the fact that T k j = h k jS k h k j,\fj / k, it follows that S k satisfies the constraints 
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given in Problem (O for the given k. Therefore, Sk must be a feasible solution for Problem (O with the given k and 

r*. 

Next, we need to prove that Sk is indeed the optimal solution of Problem © for any given k, and thus the 
corresponding achievable rate Rk is equal to the optimal value of Problem (©, which is Ck(Tk)- We prove this result 
by contradiction. Suppose that the optimal solution for Problem (O, denoted by S%, is not equal to Sk for a given k. 
Thus, we have 

R k <lo S (l + j" kS } hkk , ) (3D 



Ylj^k Fjk + °k 

}2j/l, h fk S 3 h jk + °l t 



log ( 1 + ^ "Tr. . o I = T h . (32) 



Furthermore, since h^S\hkj < Tkj,Vj ^ k, we have for any j ^ k, 



^=log 1 + - ; - (33) 



< log ( 1 + x? ; /; ^ j (35) 

V F ij + h kj S k h kj + Vj J 
= log | 1 + ^jgAj \ A (36) 

Thus, for another set of transmit covariance matrices given by Si, ... , S^-i, S\, Sk+i, • • • , S^, the corresponding 
achievable rate-tuple for the MISO-IC, (n, . . . ,rjf), satisfies that > Rk and r 3 - > Rj,Vj ^ k, which contradicts 
the fact that {R\, ■ ■ ■ , Rk) is a Pareto-optimal rate-tuple for the MISO-IC. Hence, the presumption that Sk / S* k for 
any given k cannot be true. Thus, we have Sk = S\ and Rk = Ck(Tk),Vk. Proposition 13.21 thus follows. 

Appendix C 
Proof of Proposition 14. II 

As given in Proposition 14. 1 1 with T, the corresponding optimal values of the problems in © for all fc's, Cfc(rjt)'s, 
correspond to a Pareto-optimal rate-tuple for the MISO-IC, denoted by (Ri, . . . , Rk)- Let Si, ... , Sk denote the set 
of optimal solutions for the problems in ©. We thus have 

C k (T k ) = R k = log (l + h ^ khkk ) , k = l,...,K. (37) 

V E j7 Lk T jk + ^kJ 

Next, we prove Proposition 14. 1 1 by contradiction. Suppose that there exists a pair of (i,j) with \Dij\ / 0, where 

is defined in (fT"3T ). Define a new T' over T, where all the elements in T remain unchanged except [Ty, Tji] T being 

replaced by 

[T^, T'^ = [T tjl T ji ] T + 5 tj -d l3 (38) 
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where Sij > is a small step-size, and dij is any vector that satisfies Dijdij > (component- wise), with one possible 
value for such dij is given by (fT6l ) in the main text. With T', the optimal solutions for the problems in (O remain 
unchanged \/k / while for those with k = i and k = j, the optimal solutions are changed to be S* and S*, 
respectively. Accordingly, the new achievable rates in the MISO-IC for any k ^ i,j are given by 



T k = log 11 + 

= log ^1 + 

>R k 



T,i^k,i,j hf k Sih lk + h ik S*hik + hf k Sjh jk + a 
h kk S k h kk \ 

£#M,j + h fk S t h ik + h fk S *j h jk + 



(39) 

(40) 
(41) 



where (HTb is due to d3"7T ) and the facts that hflS*h ik < T ik and h^ k S*hj k < Tj k . Also, it can be shown that 



H c* i 



log 1 + 



Ei^jhilSihu + hftSpji + ai 
.He?* 



log 1 + 



>c*(it 



Ysl ti.j r H + h^Sjhji + cr; 



(42) 

(43) 
(44) 



where (1441 is due to the facts that h^Sjhji < T'^ and S 1 * achieves the optimal value of Problem © with k = i and 
the given 1^, denoted by CiiT'j). Similarly, it can be shown that rj > Cj(T'j). Thus, from (l38l ) and Dijdij > 0, it 
follows that with sufficiently small <5y 

' CiiT'i) 
CjiT'j) 

c,(r 4 ) 

CJf(Ti) 
FL 



> 



(45) 



(46) 



> 



i? i 



(47) 



Therefore, we have found a new set of achievable rate-tuple for the MISOTC with T', (n, . . . , tk), which has 
Ti > Ri, rj > Rj, and r k > R k ,Vk ^ Clearly, this contradicts the fact that (Ri, . . . ,Rk) is Pareto-optimal for 
the MISO-IC. Thus, the presumption that there exists a pair of (i, j) with \Dij\ ^ cannot be true. Proposition 14.11 
thus follows. 
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